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Abstract 

A partial wave analysis is presented of J/ip — > <Ptt + tt~ and <pK + K~ from a sample of 
58M J/ip events in the BES II detector. The /o(980) is observed clearly in both sets 
of data, and parameters of the Flatte formula are determined accurately: M = 965 ± 
8 (stat) ±6 (syst) MeV/c 2 , gi = 165 ± 10 ± 15 MeV/c 2 , g 2 / gi = 4.21 ± 0.25 ± 0.21. 
The 6tttt data also exhibit a strong tttt peak centred at M = 1335 MeV/c 2 . It may 
be fitted with /2(1270) and a dominant + signal made from /o(1370) interfering 
with a smaller /o(1500) component. There is evidence that the /o(1370) signal is 
resonant, from interference with /2(1270). There is also a state in tttt with M = 
1790l3Q MeV/c 2 and T = 270±™ MeV/c 2 ; spin is preferred over spin 2. This 
state, /o(1790), is distinct from /o(1710). The cf>KK data contain a strong peak 
due to /2(1525). A shoulder on its upper side may be fitted by interference between 
/o(1500) and / (1710). 



PACS: 13.25.Gv, 14.40.Gx, 13.40.Hq 
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The processes J /if) — > 07r + 7r~ and <f>K + K~ have been studied previously in the 
Mark III [1] and DM2 [2] experiments. Here we report BESII data on these 
channels with much larger statistics from a sample of 58 million e + e~ — > J /if) 
interactions. The / (980), / (1370) and a state with mass at 1790 MeV/c 2 and 
with spin preferred over spin 2, called the / (1790) throughout this paper, 
are studied here. A particular feature is that /o(1790) — > mr is strong, but 
there is little or no corresponding signal for decays to KK. This behavior is 
incompatible with / (1710), which is known to decay dominantly to KK; this 
indicates the presence of two distinct states, / (1710) and / (1790). 

A detailed description of the BESII detector is given in Ref. [3]. It has a 
cylindrical geometry around the beam axis. Trajectories of charged particles 
are measured in the vertex chamber (VC) and main drift chamber (MDC); 
these are surrounded by a solenoidal magnet providing a field of 0.4T. Photons 
are detected in a Barrel Shower Counter (BSC) comprized of a sandwich array 
of lead and gas chambers. Particle identification is accomplished using time-of- 
flight (TOF) information from the TOF scintillator array located immediately 
outside the MDC and the dE/dx information from the MDC. 

Events must have four charged tracks with total charge zero. These tracks are 
required to lie well within the MDC acceptance with a polar angle 9 satisfying 
| cos#| < 0.80 and to have their point of closest approach to the beam within 
2 cm of the beam axis and within 20 cm of the centre of the interaction 
region along the beam axis. Further, events must satisfy a four-constraint 
(4C) kinematic fit with x 2 < 40. 

Kaons, pions, and protons are identified by time-of-flight, dE/dx, and also 
by kinematic fitting. The a of the TOF measurement is 180 ps. Kaons may 
be identified by TOF and dE/dx up to a momentum of 800 MeV/c. The 4C 
kinematic fit provides additional good separation between (pun and <f)KK; 
residual crosstalk between these channels is negligible. 

The K + K~ invariant mass distributions for J/ip — > K + K~7t + tt~ and J/if> — > 
K + K~ K + K~ are shown in Figs. 1(a) and (c); in the latter case, the K + K~ 
pair with invariant mass closest to the is plotted. The peaks of the lie at 
1019.7 ± 0.2 and 1020.0 ± 0.2 MeV/c 2 in (a) and (c), in reasonable agreement 
with the value of the Particle Data Group (PDG) [4]. In both cases, there is 
a clear signal over a modest background of events due to K + K~7t + tt~ or 
K + K~K + K~ without a 0. The curves in (a) and (c) show the background, 
assuming it follows a phase space dependence on M(K + K~). The resulting 
background is (19.0 ± 1.5)% in (a) and (6.2 ± 1.6)% in (c). Events containing 
a are selected by requiring at least one kaon identified by TOF or dE/dx 
and \M K+K - - M | < 15 MeV/c 2 . 

Before discussing the main physics results, it is necessary to deal with an 
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important background arising in J/ifi — > K + K~n + n~ . Events for the study of 
this background channel are selected in a sidebin having M(K + K~) = 1.045- 
1.09 GeV/c 2 . Fig. 2 shows Dalitz plots and mass projections for these sidebin 
events; Dalitz plots for (fin + n~ and (fiK + K~ data are shown in Fig. 3. For the 
K + K~ n + n~ sidebin, there is a strong peak in the (fin mass distribution of 
Fig. 2(b) centred at 1500 MeV/c 2 with a full-width of 200 MeV/c 2 . This (fin 
peak is of interest because of an earlier report of a possible exotic state close 
to this mass with quantum numbers J p = 1~ [5]. The reflection of this peak 
produces a horizontal band at the bottom of Fig. 2(a); it projects to a broad 
peak centred at 2450 MeV/c 2 in Fig. 2(b). For K + K~K + K~ sidebin events 
of Fig. 2(c), there is no corresponding peak at low mass in (fiK, Fig. 2(d). 

In order to investigate the nature of this peak, we select events in the mass 
range 1400-1600 MeV/c 2 from Fig. 2(b). Mass distributions of K+n~ and 
K~n + pairs are shown in Fig. 2(e) and corresponding distributions for K + n + 
and K~n~ in Fig. 2(f). There is a strong K*(890) peak visible in Fig. 2(e) but 
only a broad peak in Fig. 2(f). It can be shown that the presence of K*(890) 
in the background, combined with kinematic selection in a narrow range of 
K + K~ masses, can generate the peak position and width of the spurious peak 
in (fin. 

A similar effect arises in selected (finn events. Fig. 3(b) shows M((fin) for events 
selected as (fin + n~ by requiring M(K + K~) within ±15 MeV/c 2 of the (fi mass. 
There is again a (fin peak, centred now at 1460 MeV/c 2 . Again it can be shown 
that the peak is consistent entirely with background. There is no significant 
evidence for an exotic (fin state. If it were misinterpreted as a (fin state, fits 
show that it requires a (fi combined with an L = 1 pion coming from the 
7^(890), hence quantum numbers J p = 1~, 0~, or 2~. 

We have carried through a full partial wave analysis in three alternative ways: 
(a) making a cut in order to exclude events lying within ±80 MeV/c 2 of 
the central mass of i^*(890), which is slightly narrower than the selection of 
Fig. 1; (b) including into the fit an incoherent background from K*(890)Kn; 
and (c) making a background subtraction which allows for the shift in mass 
and width between sidebin and data for the background peak in (fin at 1500 
MeV/c 2 . Results of these three approaches agree within errors. We regard the 
first method as the most reliable, since it is independent of any modelling 
of the background. Figs. (4) and (5) show the fit from this approach. The 
cut against K*(890) eliminates the (fin peak at 1460 MeV/c 2 , as shown in 
Fig. 4(d). It also eliminates backgrounds due to channels JT*(1430)fr*(890), 
observed in the final state K + K~n + n~ . It reduces the background under the 
(fi in Fig. 1(e) to (13.5 ± 1.4)%; after the background subtraction, the number 
of <fin + n~ events falls to 4180. 

The branching fractions for production of (finn and (fiKK are determined allow- 
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ing for the efficiencies for detecting the two channels and correcting for unob- 
served neutral states. Results are: B(J/ip — > (fmir) = (1.63±0.03±0.20) x fCT 3 
and B(J/ip -> <pKK) = (2.14 ± 0.04 ± 0.22) x 10~ 3 . The main contributions 
to the systematic errors come from differences between data and Monte Carlo 
simulation for the selection, if* (890) cut, and particle identification; from 
uncertainties in the MDC wire resolution; and the total number of J/ip events. 

We turn now to the physics revealed by diagonal bands in the Dalitz plots 
of Fig. 3 and mass projections of Figs. 4 and 5. There is a strong / (980) — > 
7T + 7r - signal in Fig. 4(c) and a low mass peak in Fig. 5(c) due to / (980) — > 
K + K~ . Secondly, the 07r + 7r~ data exhibit in Fig. 4(c) the clearest signal yet 
observed for /o(1370) — > 7r + 7r~. Several authors have previously expressed 
doubts concerning the existence of / (1370), but present data cannot be fitted 
adequately without it. Both Mark III and DM2 groups observed a similar 
peak with lower statistics. There have been earlier reports of similar but less 
conspicuous peaks in ttti — > KK from experiments at ANL [6,7] and BNL [8]. 
A third feature in the 07r + 7r~ data in Fig. 4(c) is a clear peak at around 1775 
MeV/c 2 . 

The (f)K + K~ data of Fig. 5(c) contain a strong ^(1525) peak. However, it is 
asymmetric and may only be fitted by including on its upper side /o(1710) 
interfering with other components. 

We now describe the maximum likelihood fit to the data. Amplitudes are 
fitted to relativistic tensor expressions which are documented in Ref. [9]. The 
full angular dependence of decays of the <fi and 7r + 7r~ or K + K~ resonances 
is fitted, including correlations between them. The line-shape of the <fi is not 
fitted, because the is much narrower than the experimental resolution. We 
include production of J p = + resonances with orbital angular momentum 
£ = and 2 in the production process J/ip — > <j)f . For production of / 2 , there 
is one amplitude with £ = and three with £ — 2, where £ and the spin of 
f'2 may combine to make overall spin S — 0, 1 or 2. The one possible £ = 4 
amplitude makes a negligible contribution. The acceptance, determined from a 
Monte Carlo simulation, is included in the maximum likelihood fit. All figures 
shown here are uncorrected for acceptance, which is approximately uniform 
across Dalitz plots except for the effect of the if* (890) cut. 

The background subtraction is made by giving data positive weight in log 
likelihood and sidebin events negative weight; the sidebin events (suitably 
weighted by K + K~ phase space) then cancel background in the data sample. 

The 07r + 7r~ and <f)K + K~ data are fitted simultaneously, constraining reso- 
nance masses and widths to be the same in both sets of data. Table 1 shows 
branching fractions of each component, as well as the changes in log likelihood 
when each component is dropped from the fit and remaining components are 
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re-optimised. 



We begin the discussion with 4>K + K~ data. There is a conspicuous peak 
due to /2(1525). The shoulder on its upper side is fitted mostly by / (1710) 
interfering with / (1500), but there is also a possible small contribution from 
/o(1790) interfering with / (1500). The overall contributions to <fiK + K~ are 
shown by the upper histograms in Figs 5(c) and (d). 

The / 2 (1270) signal reported below in 07T7t data allows a calculation of the 
/2(1270) — > K + K~ signal expected in <f)K + K~, using the branching fraction 
ratio between KK and nn of the PDG. Its contribution is negligibly small. 

We discuss next the fit to / (980). In 07r + 7r~ data, it interferes with a broad 
component well fitted by the a pole [10]. This component interferes construc- 
tively with the lower side of the / (980) in Fig. 4(c). Its magnitude is shown 
by the lower curve in Fig. 4(e). 

The / (980) amplitude has been fitted to the Flatte form: 



M 2 - s- i(gip^ + 92PKR) ' 

Here p is Lorentz invariant phase space, 2k/y/s, where k refers to the 7r or K 
momentum in the rest frame of the resonance. The present data offer the op- 
portunity to determine the ratio ^2/^1 accurately. This is done by determining 
the number of events due to /o(980) — > nn and — > K + K~ and comparing with 
the prediction from the Flatte formula, as follows. After making the best fit to 
the data, the fitted /o(980) — > 7r + 7r~ signal is integrated over the mass range 
from 0.9 to 1.0 GeV/c 2 . The fitted /o(980) — > K + K~ signal is integrated over 
the mass range 1.0-1.2 GeV/c 2 , so as to avoid sensitivity to the tail of the 
/o(980) at high mass. The latter integral is given by 

0.5 J ds\f(980)\ 2 p(K + K~)e(K + K-) (2) 

and the former by 

ds\f(980)\ 2 p(nn)e(n + n-). (3) 

Here e(K + K~) and e(7r + 7r~) are detection efficiencies. The numerical factors 
at the beginning of each expression take into account (a) there are equal 
numbers of decays to K + K~ and K°K° and (b) two-thirds of irn decays are 
to 7r + 7r~ and one third to 7r°7r°. 
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By an iterative process which converges rapidly, the ratio 92/91 is adjusted 
until the ratio of these two integrals reproduces the fitted numbers of events 
for <pK + K- and <pn + it~. The result is g 2 /g 1 = 4.21 ±0.25 (stat) ±0.21 (syst). 
The systematic error arises from (i) varying the choice of side bins and the 
magnitude of the background under the peak, (ii) changes in the fit when 
small amplitudes such as / (1500) and / 2 (1270) — > K + K~ and a — > K + K~ 
are omitted from the fit, (iii) changing the mass and width of other components 
within errors and different choices of a parameterization from Ref. [10]. The 
result is a considerable improvement on earlier determinations. The mass and 
<7i are adjusted to achieve the best overall fit to the peak in 07r + 7r~ data. 
Values are M = 965 ± 8 ± 6 MeV/c 2 , g 1 = 165 ± 10 ± 15 MeV/c 2 . 

The ratio 92/91 is only weakly correlated with M and g\. However, g 2 is rather 
strongly correlated with M. This arises because the term ig2PKK{s) in the 
Breit-Wigner denominator, eqn. (1), continues analytically below the KK 

threshold to —g 2 \J (AMj^/ s) — 1. It then contributes to the real part of the 
Breit-Wigner amplitude and interacts with the term (M 2 — s). We find that 
the correlation is given by dg2/dM = —5.9; the mass goes down as g 2 goes up. 
Other correlations are weak: dgi/dM = —0.75 and dr/dgi = —0.068, where 
r = 92/91- 

We consider next the peak in 07T7T centred at a mass of 1335 MeV/c 2 . An 
initial fit was made to /2(1270) and one fa. The /o optimizes at M — 1410 ± 
50 MeV/c 2 , T = 270 ± 45 MeV/c 2 , where errors cover systematic variations 
when small ingredients in the fit are changed. However, both / (1500) and 
/o(1370) can contribute. Adding / (1500), log likelihood improves by 51: an 
8.5 standard deviation improvement for four degrees of freedom. Also the fit 
to the 7T7T mass distribution improves visibly. Therefore three components are 
required in the 1335 MeV/c 2 peak: / 2 (1270), / (1370) and / (1500). Removing 
/o(1370) makes log likelihood worse by 83, a 10.8 standard deviation effect. 

Angular correlations between decays of <fi and / 2 are very sensitive to the pres- 
ence of / 2 (1270), which is accurately determined. It optimizes at M — 1275±15 
MeV/c 2 , T = 190 ±20 MeV/c 2 , values consistent with / 2 (1270). The fact that 
its mass and width agree well with PDG values rules out the possibility that 
the remainder of the signal in this mass range is due to spin 2; otherwise the fit 
to / 2 (1270) would be severely affected. Angular distributions for the remaining 
components are indeed consistent with isotropic decay angular distributions 
from spin 0. 

The / (1370) interferes with / (1500) and / 2 (1270). This helps to make / (1370) 
more conspicuous than in other data. However, because of the interferences, 
its mass and width are not accurately determined. The mass of /o(1370) is 
1350 ± 50 MeV/c, where the error is the quadratic sum of the systematic and 
statistical errors. 
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The width of /o(1370) is somewhat more stable. It is determined essentially by 
the full width of the peak in 07r + 7r~ of 270 MeV/c 2 ; interferences with / 2 (1270) 
and / (1500) affect this number only by small amounts and the fitted width 
is 265 ±40 MeV/c 2 . If both / (1370) and / (1500) are removed, log likelihood 
is worse by 595. Removing / (1500) from the fit perturbs the mass fitted to 
/o(1370) upwards to 1410 ± 50 MeV/c 2 ; this is obviously due to the fact that 
/o(1370) is trying to simulate the missing / (1500) component. 

The presence of a peak due to / (1370) is strongly suggestive of a resonance. In 
order to check for resonant phase variation, we have tried replacing the ampli- 
tude by its modulus, without any phase variation. In this case, log likelihood 
is worse by 39, nearly a 9 standard deviation effect for a change of one degree 
of freedom. The conclusion is that the /o(1370) peak is resonant. It is not 
possible to display the phase directly, since it is determined by interferences 
between two / (1370) and four / 2 (1270) amplitudes. 

The magnitude of the signal due to /o(1370) — > K + K~ in the fit gives a 
branching fraction ratio 

fg|lg = 0.08 ±0.08. (4) 
B[/ (1370) -. tttt] ^ ' 



This value is somewhat lower than reported by the Particle Data Group [4]. 
The reason is the conspicuous signal in nn but absence of any corresponding 
peak in K + K~ . 

Next we consider the peak in 7r + 7r~ at 1775 MeV/c 2 in Fig. 4(c). It fits well 
with J p = 0+ with M = 17901^ MeV/c 2 , T = 270l|g MeV/c 2 . The fitted 
mass is in reasonable accord with the / (1770) reported in Crystal Barrel data 
on pp _^ ( w ) 7r o [ii] : M — 1770 ±12 MeV/c 2 , T = 220 ±40 MeV/c 2 . Allowing 
for the number of fitted parameters, / (1790) is more than a 15cr signal. It 
cannot arise from / (1710), since the magnitude of /o(1710) — > K + K~ is small 
(see Table 1), and it is known that the branching fraction ratio of /o(1710) 
between nn and KK is < 0.11 at the 95% confidence level [12]; accordingly, 
the / (1710) — > 7T7T signal in present data should be negligibly small. 

We now consider possible fits with an / 2 instead. The decay angular distri- 
bution in this mass range is consistent with isotropy. So there is no positive 
evidence for spin 2. However, four spin 2 amplitudes are capable of simulating 
a flat angular distribution. In consequence, spin 2 gives a log likelihood which 
is worse than spin by only 4.5 after re-optimising its mass and width. If 
/o(1710) is then added with PDG mass and width, it improves log likelihood 
by a further 2.0; this confirms the result from uK + K~ data that /o(1710) has 
a negligible decay to nn. Our experience elsewhere is that using four helicity 
amplitudes instead of 2 adds considerable flexibility to the fit. The spin 2 am- 
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plitude with £ = has a distinctive term 3 cos 2 a n — 1, where is the decay 
angle of the n + in the resonance rest frame, with respect to the direction of the 
recoil 0. Simulation of spin requires large J = 2, £ = 2 and 4 amplitudes to 
produce compensating terms in sin 2 a n and hence a flat angular distribution. 
Large contributions from £ = 2 are unlikely in view of the low momentum 
available to the resonance and the consequent £ = 2 centrifugal barrier. If the 
J = 2 hypothesis is fitted only with £ — 0, log likelihood is worse by 95 than 
for spin 0. We conclude that the state is most likely spin zero. 

It is not possible to fit the shoulder in <J)K + K~ at 1650 MeV/c 2 accurately by 
interference between / (1500) and /o(1790), using the /o(1790) mass and width 
found in (pirn data. Even if one accepts the poor fit this gives, the branching 
fraction ratio KK /titx assuming only one fo resonance here is 0.55 ±0.10. This 
is a factor 14 lower than that reported in Ref. [12] for / (1710). For a resonance, 
branching fractions must be independent of production mechanism. The large 
discrepancy in branching fractions implies the existence of two distinct states 
at 1710 and 1790 MeV/c 2 , the former decaying dominantly to KK and the 
latter dominantly to nn. The / (1790) is a natural candidate for the radial 
excitation of / (1370). There is earlier evidence for it decaying to An in J/ip — > 
7(47r) data [13,14], with mass and width close to those observed here. There, 
spin was preferred strongly over spin 2. 

The shoulder in <fiK + K~ at 1650 MeV/c 2 is fitted with interference between 
/o(1500) and / (1710), which is known to decay strongly to KK. If both 
/o(1710) and / (1790) are included in the fit, there is only a small improvement 
from / (1790). 

Masses, widths and branching fractions are given in Table I. The errors arise 
mainly from (i) varying the choice of side bins and the magnitude of the back- 
ground under the peak, (ii) adding or removing small components such as 
/o(1500), / 2 (1270) -> K+K~, and a -> K + K~ and (iii) varying the mass and 
width of every component within errors and using different a parameteriza- 
tions reported in Ref. [10]. It also includes the uncertainty in the number of 
J ftp events and the difference between two alternative choices of MDC wire 
resolution simulation. 

Finally, angular distributions for both production and decay have been exam- 
ined for each resonance peak. There are no significant discrepancies between 
data and fit. A fit is shown for the / (1790) peak in Fig. 6 to the decay angle 
a n of the rnr pair, with respect to the recoil 0; the deep dip at cosa^ = ±0.75 
is due to the if* (890) cut. The remaining angular distribution fits well to spin 
0. 

It is remarkable that (pmr data contain large signals due to several states which 
are predominantly non-strange: / 2 (1270), / (1370), / (1500) and / (1790); di- 
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1790_ 30 
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1521 ±5 


77 ± 15 




7.3 ± 1.1 


440 


/o(1710) 


PDG 


PDG 




2.0 ±0.7 


64 



Table 1 



Parameters of fitted resonances and branching fractions for each channel; improve- 
ments in S = log likelihood when the channel is added. PDG means that the mass 
and width are fixed to the PDG value. For the /o(980), see the parameterization in 
the text. The errors are the statistical and systematic errors added in quadrature. 

rect production with the should favour ss states. There is no agreed expla- 
nation. 

In summary, the data reported here have three important features. Firstly, the 
parameters of / (980) are all well determined. Secondly, there is the clearest 
signal to date of /o(1370) — > 7r + n~; a resonant phase variation is required, 
from interference with / 2 (1270). Thirdly, there is a clear peak in at 1775 
MeV/c 2 , consistent with / (1790); spin 2 is less likely than spin 0. If the 
/ (1790) resonance is used to fit the shoulder at 1650 MeV/c 2 in <pK + K~, the 
branching fraction to pions divided by that to kaons is inconsistent with the 
upper limit for the ratio observed in Ref. [12] for /o(1710), this requires two 
distinct resonances /o(1790) and / (1710). 
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Mass(K + K~) Mass(I?K) 




Mass(K + K) Mass(n + n) 

Fig. 1. The K + K~ invariant mass distributions for (a) J/tp — > K + K~7r + ir~ , (c) 
J/ip — > K + K~ K + K~ '; curves show the fitted background and a Gaussian fit to the 
4>; (b) and (d) show mass projections for events selected within ±15 MeV/c 2 of the 
4>; (e) and (f) show mass projections after cutting events within ±100 MeV/c 2 of the 
central mass of K*(890); curves in (e) show the fitted background and a Gaussian 
fit to the (p. 
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Fig. 2. Dalitz plots for a sidebin to the <fi with M(K + K~) = 1.045 to 1.09 GeV/c 2 for 
(a) K + K~tt + tt~ data, (c) K + K~~K + K~ data; (b) and (d) show projections against 
(fin and (fiK mass; the mass distributions for (e) M(K + n~) and M(K"tt + ), (f) 
M{K +/ n + ) and M(K~ir~) for events from (a) in the (fin mass range 1400-1600 
MeV/c 2 . 
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Fig. 3. Dalitz plots for (a) J/tp -> ^vr+vr", (c) J/i/) -» 4>K + K-; (b) and (d) show 
the projections against 07r and <j)K mass. 
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Fig. 4. (a) and (b) show measured and fitted Dalitz plots for J ftp — > <j)7r + ir~ after 
cutting i<C*(890) events, (c) and (d) show mass projections; the upper histogram 
shows the maximum likelihood fit and the lower one shows background; (e) shows the 
fo contribution from the fit (full histogram) and the lower curve the a contribution, 
(f) the / 2 (1270) contribution. 



15 




Fig. 5. (a) and (b) show measured and fitted Dalitz plots for J/ip — > <pK + K~. 
(c) and (d) Mass projections for J/ip — > <pK + K~ data, compared with histograms 
from the fit; (e) shows the /o contribution from the fit (full histogram) and the 
lower curve the /o(980) contribution, (f) the /2(1525) contribution from the fit (full 
histogram) . 
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Fig. 6. Angular distributions in </>7r + 7r~ data (crosses) for a n , the angle of the ir + 
from fj decay with respect to the direction of fj in its rest frame. The upper 
histograms shows the fit, and the lower one the background. The dashed histogram 
shows the acceptance. 
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